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Abstract 

The evolution of linear cosmological perturbations in modified theories of gravity is 
investigated assuming the Palatini formalism. It has been discussed about the stability 
problem in this model based on the equivalence between f(R) gravity and the scalar tensor 
theory. However, we study this problem in the physical frame where the matter is minimally 
coupled. In general, the stability of the superhorizon metric evolution depends on models. 
We show that the deviation from the superhorizon metric evolution is null for a specific 
choice for the nonlinear Einstein-Hilbert action, f{R) ~ -R", where n 7^ 0, 2, 3. Thus the 
stability of metric fluctuation is guaranteed in these models. We also study the matter 
density fluctuation in the general gauge and show the differential equations in super and 
sub-horizon scales. 



1 Introduction 



The discovery of the present accelerated expansion of the Universe [1, 2] can be accounted 
for either by the existence of a homogeneous component of energy with a negative pressure, 
dubbed dark energy [3], or by a modification of gravitational action with a general function 
of the scalar curvature instead of the standard Einstein-Hilbert term, named f(R) gravity 
[4]. 

Alternative field equations depend not only on the choice of action but also on the 
variational principles [5]. The Palatini formalism where the metric and the connections 
are treated as independent variables and the energy momentum tensor does not depend 
on the independent connection leads to a different theory from what is obtained from the 
metric formalism. The Palatini formalism of f(R) gravity results in second order differential 
equations due to the algebraic relation between the curvature scalar and the trace of the 
energy momentum tensor. 

Cosmological background solutions have been studied for various gravitational La- 

grangian in the Palatini formalism [6]. Their validity as cosmological theories have been 
tested for observation [7]. Also the cosmological perturbations have been investigated [8]. 

It has been known that the Palatini formalism to f(R) gravity is stable for the curva- 
ture scalar perturbation [9]. The Newtonian limit of models using the Palatini variational 
principle gives contradicting results [10]. However, the approach to this problem is using 
the equivalence between f(R) gravity and scalar-tensor theory. We will use the evolution 
of linear perturbations in f(R) models in the physical frame as already done in the metric 
formalism [11]. 

We study the evolution of linear perturbations in Palatini f(R) gravity in the physical 
frame where the matter is minimally coupled. We consider the stability of metric fluctua- 
tions at high curvature to see the agreement with high redshift cosmological observations. 
We also investigate the evolution of the matter density fluctuation. 

In the next section we review the Palatini f(R) gravity. We investigate the hnear per- 
turbation of f(R) gravity in the section III. Compared with previous works [8], we do not 
specify the gauge to find the matter density fluctuation. In section IV, we derive the sta- 
bility equation of metric fluctuations in the high curvature limit and show the stability in 
a specific model. We show the evolutions of the metric fluctuation and the density contrast 
in the superhorizon and the subhorizon scales in section V. We reach our conclusions in 
section VI. We also provide the detail calculation in the appendix. 
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2 Palatini f(R) gravity 



We consider a modification to the Einstein-Hilbert action assuming the Palatini formahsm, 
where the metric g^i, and the torsionless connection F^^ are independent quantities and the 
matter action depends only on metric 



S 



df^x^f—g 



(2.1) 



where are matter fields. Then the Ricci tensor is defined solely by the connection 

p "pa -pa _i_ pa p/3 pa p/3 (n 

whereas the scalar curvature is given by 

R = g^^^R^, . (2.3) 



We can derive the field equation of f(R) gravity in the Palatini formalism from the above 
action (2.1) 

F{R)R,, - ^g^J{R) = K^T^. , (2.4) 
where F{R) — df{R)/dR and the matter energy momentum tensor is given as usual form 



2 5{y^Cm) 



Prom the above equations we can get the generalized Einstein equation 



(2.5) 



G,. = G,. + ^-Iv^FV.F - ^V^V.F + g,^^DF - ^g^^l^{dFf , (2.6) 

where (7^,^ = R^i, — \g^vR and G^i, = R^j^,, — ^g/^uR- For the later use, we reexpress the 
above equation again as 

FG^, = - ^jV^FV^F + V^V.F + ^(/ - FR)g^, + ^g^.aF . (2.7) 

Prom the above equation we can derive 00-component of the modified Einstein equation 



F'^ 2 



F" / H\F' 

+ + 



2F ' 



(2.8) 



where primes denote derivatives with respect to In a. We also have the following useful 
equation 



-2 



H' 



K 



H FH^ 



{P + P) + 



F" /_ H\F' _3F''' 



(2.9) 
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3 Linear perturbation in f(R) gravity (In Conformal 
Newtonian Gauge) 



First we start from the metric in the conformal Newtonian gauge. The hne element is given 

by 



(l + 2*(t, f))dT^ + (l - 2^(t, x))dx'dxi 



(3.1) 



The main modifications for viable models with stable high curvature limits happen well dur- 
ing the matter dominated epoch and we can take the components of the energy momentum 
tensor as 

To° = -p(l + <5), Tf = pd,q, T; = 0. (3.2) 
From the previous equation (2.7), we can find the perturbed Einstein equation 



2 F 



4 F 



2 F2 



(3.3) 



where we use the equation (7.7) and Sf{R) — F{R)SR. If we consider the ij-component of 
the perturbed equation, then we can find 



$ - * = — , 
F ' 



(3.4) 



where we assume the null anisotropic stress. If we use the above equation (3.4), then we 
can express the other components of the perturbed Einstein equation (3.3) 



3H' 



, IF' , „ AF" IH'F' IF' H' x 
$' + ^' + + ( \ \ \ ^ 1]$ 

^2F^ ^ V2F 2F^ 2HF 2F H ) 



IF" 



711 



/ r- IH'F' 3F' H' s 

+ ( \ \ hi)* 

V2F F2 2H F 2F H ^ 



H 



1 F' 

$' + *' + $ + * + - — ($ + *) 

2 r 



(3.5) 
(3.6) 



*" + *"+(4 + f)*'+(3^ + 4 + fK+(-^-(2 + |)5 



+7r + 3)*+(3-F + (6 + 3-)- + 3- + 3)l. 



0. 



(3.7) 



To capture the metric evolution, let us introduce two parameters as in the reference [11]: 
9 the deviation from ( conservation and e the deviation from the superhorizon metric 
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evolution 



(' + - H'q = 



+ ^' _ + ii 

//' \H H' 

where we define the dimensionless quantity 



I* 



aH 



Be, 



F H' 

Prom the above equations, we can find the expression for 6 and e, 
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(3.8) 
(3.9) 

(3.10) 



H \aH 



Be 



B' 3 H' H" 
B 2 H H' 



($-^) + 



1 



2H' 



r n 
H 



H 



H ) 



1 H'^ r.2 



Hq 



(3.11) 



aH 



Be 



H' „ H" H' 
— B + — + — + 3 
H H' H 



5 H' H" 

+ B - 6 — 

2 H H' 



H 



//2 



'B^ 



+ (4- 



H' 



H'2 

B' 

'b 



-2 



^2 



li/' 



/2 



B ^ 2H^ 



B' 



H' fH" H' N 1 

3— + 3-2(— + — + 3)- 
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H' 



H 



B 



1 H'^ r.2 . 

+2— 



B'+{2^ + 



H 



2 i/^ 



B''' f H" H' H', H'^. ron on 

V + + ¥ + + +(2^+2^ 

}i\ H" H"^ H" H' fH" H' xT 



($-*) + 
i7' 



//2 
^72 



-H^^2m 



H' H 

2H) 



2 



'fh 



Q 



(3.12) 



Prom equation (3.11), we can recover the conservation of Newtonian gauge when $ 
and F is a constant. 

In addition to these equations, we can find very useful equation from the structure 
equation (2.4). By taking the trace of this equation (2.4) and differentiate with In a, we 
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have 



R' 



F,rR - F 



-r. 



(3.13) 



Also by taking the perturbation of the equation we have 



SF = Fp5R 



From these two equations, we can find 



,R 



F.R - F 



-k'ST. 



(3.14) 



5F 1 F' 

_ = ---i = $-*. (3.15) 
Prom this equation we can find the evolution equation of matter density fiuctuation 



B' H' 



H" H' 



H' 



B 



12 



5"+ 2— + — B + 2— - — + 3U'+ \ —B'-2—- 



B H 



12 jj// jj/2 jjf 



H 



H" H' 



H' H 

B + 

3*'. 



H 



S2 



H" H' 



H" H" H"^ H'^ H' 



'b 



(3.16) 



Compared with previous works [8], we do not specify the gauge of matter density to solve 
the matter density fiuctuation. 



4 Stability of metric fluctuations 



Unstable metric fiuctuations can create order unity effects that invalidate the background 
expansion history. We can derive the evolution equation of the deviation parameter. If we 
differentiate the equation (3.12) and consider the evolution in the superhorizon scale, then 
we have 



„ / B' H' H" 

e" + 2 h — B + — 

\ B H H' 



+ 



IH' 



12 



2 



B^ 



-4} 



B 



r H" H'^ 

2 4— 

L H' 

F{^,^,Hq), 



H' 

3H' 
' 2ll 



B + Q' 



H' , 



H" H' iB' 

5 h4 h9 — 

H' H i B 



-2- 



H" 

IF 



H' 
ll 



g + 4 



H 



12 



H' 



(4.1) 



where wc use equations (3.7) and (3.9) and F(^, $) is the source function for the deviation 
e and define Q as 

H" H' 

Q = ^ + ^ + 3. (4.2) 
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The above equation is different from the metric formahsm [11]. The stabihty of e depends 
on the sign of the coefficient of the term proportional to e. In the metric formahsm e is 
stable as long as B > 0. However, the stability is complicate and need to be checked for 
each model in the Palatini formalism. 



4.1 A particular example : f{R) = /3R^ 

We demonstrate the general consideration of the previous subsection with a specific choice 
for the nonlinear Lagrangian, f{R) = where n ^ 0,2,3. The background is simply 

described by a constant effective equation of state in this model. The Hubble parameter 
scales as H'^ ~ a"^/". Then it is easy to write it with its derivatives in terms of In a 

£!=f-±y. (4.3) 

H 2n H \ 2n) ^ ^ 

Here the scalar curvature is i? = 3(3 — n)H'^/{2n). From this fact, we can also find the 
derivatives of F with respect to In a 

F' _ F" _ 3(1 - n) F" _ F'" _ ( ?,{l-n) V 

F ~ F' ~ n ' F ~ F' ~\ n ) ' ^ ' ^ 

If we use above equations (4.3) and (4.4) into (3.12), then we find that the deviation from 
the super horizon metric evolution is null, e = 0. 



5 Metric and matter density evolutions 
5.1 Super horizon evolution 

We consider the metric evolution in superhorizon sized, k/{aH) <^ 1. In this case, the 
anisotropy relation of the equation (3.11) becomes 

$ - * ~ (b + (5.1) 

where A is given by 



5(25f + 5) 

Prom the above equation (5.1), we can find the superhorizon evolution equation of $ from 
the equation (3.9) 

:^ + 2— 5 + — - — + 4-C U'+ — + — 5+ — + 3-CU~0, (5.3) 
B H H' H ) \B H H' ) '^^ 
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where C is defined as 



B' H' H" H' 

h2 — B + 2 h3 

B H H' H 



(5.4) 



B + A + 1 

If we use this equation (7.12), then wc have the evolution equation of matter density (3.16) 
B \ ^„ ( B' H' „ B'A - BA' H' B 



1 + 



B + A 



+ 1 







^ B' 























H' B'A - BA' 



B' H" 



H B + A 

w 



H" H' \ 
2^-^ + 3)5' 



H'^ iy'i 



H 
B 



B"A-BA" (B'A- BA') {B' + A') H'B'A-BA' ( H'\' B 



{B + Ay 



- 2 



{B + Ay (B + A) H {B + Ay 



+ 



H" H" H"^ H'^ H' 



-4 



H" H' 



H J B + A 
0. 



(5.5) 



5.2 Superhorizon evolution in a particular example 

Now we can check the evolution equations in the previous sTibsection in a particular case, 
f{R) ~ R^- In this case, we can simplify the following quantities 



5 = 2(n-l), A = -4(n- 1) = -2S, C = ^ = 
Prom this, we can also simplify the evolution equations (5.3) and (5.5) 



^„ 9-4n^, 



2n 



S' 



0, 
0. 



(5.6) 

(5.7) 
(5.8) 



The evolution equation of the Newtonian potential has no terms proportional to thus 
$ = constant is a solution to the equation. Also the matter density fiuctuation has the 
same for as general relativity, 5 — constant. 



5.3 Subhorizon evolution 

For subhorizon scales where k/aH ^ l,we can find the Poisson equation from the equation 
(3-5) 

A;2($ + ^)c._^!!^5. (5.9) 
F 

If we use equations (3.15) and (5.9), then we have 

/ o a^H^ F'\S F'5 
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Prom this equation we can find 

$ ~ — ^ . 

We can differentiate (7.10) and use the above equation (5.11) to get 



(5.11) 



— + B + h — + 6$'- 

B 2 H H' H 



B^ + 



H" 



H" 21 H' 



2 4 

H H' 2 H 



B + 



H 



52 

//2 



+ 



^2^ - 1 

H' 



2 — B' 

B H 



H" 



H' 



+ 



H" H' 



(5.12) 



If we differentiate the equation (5.10) and put into the equation (3.16), then we have 



( B' 3H' H" H' fSH' , S'2 



12 



^ R2 I 



3H" 3H'^ rH'lr. 



2H 

H" , H" 

-3 4 — 

H H' 



52 



H" H' ,B' 



B 



2 \ 



4+" 



H' 



B 



15 = 0. 



(5.13) 



5.4 Subhorizon evolution in a particular example 

We can use the previous relation (5.6) into the evolution equations (5.12) and (5.13) 

(5.14) 
(5.15) 

Even though the superhorizon scale evolutions of $ and 5 are same to those of general 
relativity, the subhorizon scale evolutions of them show different behaviors from those of 
general relativity as expected [12] . 



^„ ^ 3(3 -n) ^, ^ 3(14n2 + 19n-36) ^ ^ 
2n 4^2 

^,^3(2-n)^, = 
2n 



6 Conclusions 

We have analyzed the cosmological evolution of linear perturbations in Palatini f (R) gravity 
to see the stability of metric fluctuations. We have also considered the matter density 
fluctuation in the Newtonian gauge. 

Compared with metric f(R) gravity, we have shown that the stability of metric fluc- 
tuations in the high redshift limit of high curvature is not simply expressed. We need to 
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check each model for the stabihty. However, we have found that the deviation from the 
superhorizon metric evolution is null for a specific choice of the nonlinear Einstein-Hilbert 
action, f{R) ~ and stability of this model is guaranteed. 

We have investigated the evolution equations of Newtonian potential and matter density 
contrast in super and sub-horizon scales. In the specific model, superhorizon evolutions of 
Newtonian potential and matter density fluctuation are same to those of general relativity. 
However, subhorizon evolutions show the different behaviors from the general relativity 
case. This will give us the method to probe the possibility of any modification of gravity. 



7 Appendix 



From the equation (2.2) we can derive the Ricci tensor and the scalar curvature by using 
the metric relation 

R,. = R,. + l^V,FV,F - ^W^V.F -^g,,jDF , (7.1) 
R = R-3^DF + ^j-^{dFr. (7.2) 
We can rewrite field equation (2.6) as the form of Einstein equation plus corrections 

C^. = + (1 - F)R^, - ^iv^FV.F + V^V.F + ^(/ - R)g^.u + \9,..^F . (7.3) 
We can check that the correction terms are covariantly conserved by using the useful relation 



Also by taking the trace of the equation (2.4) we can find 



F" / H\F' 



+ J ■ (7.4) 



F y H ^ F 2 F2 

By differentiation of the equation (2.9) we can derive another useful equation 

F'" F"F' ( H' \F" 3F'^ ( H' 9\ F'^ ( H" H'^ H' \ F' 
-F-'^+['H+'h^-2l^-['H^-2)l^^[ir + lP^'H-')F 

The perturbed equations for the Ricci tensor and the scalar curvature are obtained from 
the equations (2.2) and (2.3) 

3 V^FV F 1 

5R^ = + ^^(V'^FV.F) - 3^^5F - -^(V'^V.F) 
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5R 



fVV^F) / 1 UF \ 

5R - piBF) + 3^SF + ^,S{dFf - 3^5F . 



(7.6) 
(7.7) 



If we differentiate tfie equation (3.6) and use tfie equation (3.7), then we have 



^f_($'_^') + 
2 ' 



3F" /9 3H\F' 



5F" /7 5H\F' ^H' 



1 I 

K p q 



(7.8) 



If we use the equation (2.9) and adopt Hq' — — ^, then we can rewrite the above equation 
(7.8) as 



B' H' „ H" 
— + — B + — + 3 
B H H' 



We can find $' from the equation (3.6) and (7.9) 

1 



TTl 

($-*) s* = o. 

H 



B' 3H' H" ^ 



1 K^p 

2FH 



(7.9) 



(7.10) 



Prom the equations (3.15) and (5.1) we can derive 



3* 



B 



B + A 



B'A - BA' H' B 



[B + Af H{B + A) 



If we differentiate this equation (7.11), then we have 



3*' 



B 



B + A 



B'A - BA' H' B 



[B + Ay 



H B + A 



S' + 



{B'A - BA') {B' + A') H' B'A - BA' 



B"A - BA" 
(B + Af 
fH'V B 



(B + Ay (B + A) H (B + Ay \H ) B + A 



(7.11) 



5. (7.12) 
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